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ABSTRACT 


A  series  of  experimental  measurements  of  the  visible  light 
angular  scattering  profiles  of  clouds  is  described.  The  results 
for  single  layer  overcast  clouds  were  found  to  correlate  well  with 
standard  hemispheric  and  narrow-angle  pyrheliometric  transmittances. 
The  resulting  expression  for  the  angular  scattering  profile  of  an  over¬ 
cast  cloud  layer  is: 


where 


T(0)  =  A  +  B  sin2G 


A  =  -0.  0056  +  0.  89  X  T 


B  =  -0.028  +  1.09  X  T 


w 


and  0  is  the  half  angle  of  the  field  of  view,  Tn  the  transmittance 
measured  with  a  narrow-angle  pyrheliometer ,  and  T^  the  transmit 
tance  measured  with  a  hemispheric  pyrheliometer. 
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I.  INTRODUCTION  AND  SUMMARY 


The  light  scattering  properties  of  clouds  are  of  great  interest  for 
a  variety  of  reasons.  In  particular,  they  have  a  major  impact  on  the 
performance  of  several  classes  of  optical  systems  for  which  clouds  may 
serve  as  a  background  or  foreground.  It  is  important  in  the  understand¬ 
ing  of  the  performance  of  these  systems  that  the  cloud  scattering  charac¬ 
teristics  be  known.  Many  experimental  and  theoretical  studies  of  cloud 
scattering  have  been  and  are  being  carried  out.  However,  one  important 
aspect,  the  angular  distribution  of  the  scattered  sunlight  that  penetrates 
the  cloud,  has  been  neglected.  In  this  paper  the  results  of  an  experimental 
measurement  of  this  distribution  are  presented.  An  analysis  of  these 
measurements  resulted  in  an  empirical  equation  that  permits  prediction 
of  this  angular  distribution  for  a  variety  of  clouds  from  standard 
pyrheliometric  measurements  made  at  most  major  weather  stations. 

In  an  earlier  study  (Ref.  1)  a  large  body  of  actinometric  data  was 
analyzed.  Probability  distributions  of  the  transmittance  of  clouds  were 
developed.  These  actinometric  data  were  measured  with  two  types  of 
pyrheliometers:  a  wide  -angle  instrument  that  accepts  radiation  from 
the  entire  celestial  hemisphere,  and  a  narrow-angle  instrument  that 
tracks  the  sun  and  accepts  radiation  from  a  5.  /°  full-angle  heliocentric 
cone.  The  wide-angle  and  narrow-angle  transmittances  were  calculated 
from  these  two  measurements  and  from  the  known  incident  solar  irradiance. 
The  value  obtained  from  the  wide-angle  pyrheliometer  is  the  effective  trans¬ 
mittance  for  a  system  that  accepts  radiation  from  a  wide  range  of  angles. 
This  transmittance  is  also  useful  in  estimating  the  light  available  for 
ordinary  photography.  On  the  other  hand,  the  transmittance  value  obtained 
from  the  narrow -angle  pyrheliometer  applies  to  a  system  that  has  a 
suitable  narrow  acceptance  angle.  The  suitability  of  the  angle  depends, 
unfortunately,  on  a  number  of  parameters,  including  the  cloud  height 
and  the  position  of  the  source  of  light  being  measured.  If  the  effective 
narrow-angle  transmittance  is  needed  for  a  range  of  conditions,  it  must 
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be  measured  as  a  function  of  the  scattering  angle.  The  experiment  des¬ 
cribed  in  this  paper  was  designed  to  perform  these  measurements. 

An  instrument  was  developed  to  measure  the  total  irradiance 
from  the  sun  am’  the  circumsolar  sky  for  a  series  of  angular  fields  of 
view  irom  1°  to  32°  haif  angle.  These  measurements  were  made  and 
recorded  automatically  every  ten  minutes  from  sunrise  to  sunset  over 
a  six  month  period, 

Eglin  Air  Force  Base  was  selected  as  the  site  for  the  experiment 
for  two  reasons.  The  first  was  that  this  geographical  area  experiences 
the  frequent  occurrence  of  suitable  cloud  formations.  The  second  was 
that  the  well-equipped  Air  Force  weather  station  at  the  site  could  provide 
concurrent  standard  meteorological  observational  data  including  radar 
cloud  height  measurements. 

A  total  of  329l  sets  of  measurements  were  made  during  the  course 
of  the  experiment.  Cf  these,  334  sets  were  obtained  when  there  was  a 
single  overcast  layer  of  clouds.  Because  of  the  difficulty  in  modeling 
mathematically  all  but  this  single-layer  overcast  data,  the  analysis 
presented  here  is  based  on  these  334  data  sets. 

A  modei  equation  that  gives  reasonable  agreement  to  the 
measured  data  is: 

T  =  A  +  B  sin“6 

where  T  is  the  transmittance,  0  the  half  angle  of  the  field  of  view,  and 
A  and  B  parameters.  These  parameters  were  found  to  correlate  well 
w'ith  the  standard  wide-  and  narrow-angle  pyrheliometric  : ’-ansmittance 
measurements.  These  parameters,  expressed  in  terms  of  Tw  and  Tn> 
which  are  the  transmittances  measured  with  the  wide-angle  and  narrow- 
angle  pyrheliometers,  respectively,  are: 
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A  =  -0.  0056  +  0.  89  X  T 

n 

B  =  -0.  028  +  1.  09  X  T 

w 

This  model  represents  an  angular  scattering  distribution  consisting  of  an 
undeviated  component  plus  a  component  which  has  been  diffusely  scattered. 
This  should  serve  as  an  adequate  model  for  cloud  scattering  profiles. 


. . .  *...  .  *x;  -  '  > 


II.  BACKGROUND  AND  PREVIOUS  WORK 


The  optical  engineer  needs  a  fairly  detailed  knowledge  of  the  transmission 
and  scattering  properties  of  clouds  when  designing  a  visible  optical  system  to 
operate  in  the  presence  of  clouds.  While  this  subject  has  been  treated  exten¬ 
sively  from  a  theoretical  view  point  (Ref.  2),  surprisingly  few  experimental 
measurements  have  been  made  by  which  the  theories  can  be  checked. 

Haurwitz  (Ref.  3)  has  correlated  measurements  of  the  total  radiation 
from  the  sun  and  sky  with  cloud  type  and  solar  elevation  angle.  Unfortunately, 
these  measurements  give  no  information  about  the  angular  scattering  distribu¬ 
tions  produced  by  the  clouds  since  the  measurements  were  made  with  a  2n 
field  of  view  detector.  Gibbons  and  coworkers  (Refs.  4-6)  have  made  studies 
of  the  effect  of  field  of  view  on  horizontal  atmospheric  transmission  measure¬ 
ments.  These  measurements  give  the  scattering  properties  of  the  atmosphere 
plus  aerosols,  haze,  etc.,  but  tell  nothing  about  clouds.  Arnulf  and  coworkers 
(Ref.  7)  have  made  transmission  measurements  of  haze  and  fog  but  unfortunately 
only  with  a  detector  with  a  narrow  angle  field  of  view.  The  aircraft-based 
measurements  of  Neiburger  (Ref.  8)  and  Griggs  and  Marggraf  (Refs.  9  and  10) 
of  the  transmission  of  solar  radiation  through  clouds  were  again,  unfortunately  , 
made  only  with  2  it  field  of  view  detectors. 
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III.  INSTRUMENTATION  AND  EXPERIMENTAL  DESIGN 


A.  EXPERIMENTAL  DESIGN 

In  designing  the  experiment,  several  requirements  were  laid  down. 

The  first  was  that  the  experiment  should  run  automatically  with  a  minimum 
of  maintenance  and  adjustment.  This  was  necessary  because  a  large  number 
of  measurements  would  need  to  be  made  in  order  to  accumulate  a  good 
statistical  sample. 

The  second  requirement  was  that  the  experimental  site  be  subject 
to  a  wide  variety  of  cloud  conditions.  The  panhandle  region  of  Florida 
experiences  the  tail  ends  of  the  winter  storms  that  sweep  down  across  the 
United  States  from  Canada  and  the  Arctic.  In  the  summer,  cumulus  clouds 
are  prevalent  due  to  the  moisture -laden  Gulf  air. 

The  third  requirement  was  that  the  experimental  site  be  located  near 
a  well-equipped  weather  station.  This  would  facilitate  concurrent  observations 
of  meteorological  conditions  on  a  routine  basis. 

The  fourth  requirement  was  for  simultaneous  narrow-  and  wide-angle 
pyrheliometric  measurements.  This  would  allow  for  testing  possible  correla¬ 
tions  of  the  expe -imental  measurements  with  standard  solar  radiation  measure¬ 
ments  that  are  made  on  a  routine  basis  at  many  stations  around  the  world. 

The  first  requirement  was  met  in  the  actual  design  of  the  instrument. 
The  second  and  third  were  met  in  the  experimental  site  selection.  Eglin  Air 
Force  Base  near  Pensacola  in  no rthwe stern  Florida  (latitude  30°29!N,  longi¬ 
tude  86°31'W)  has  a  well-equipped  weather  station  staffed  by  personnel  of 
Detachment  10  of  the  6th  Weather  Wing,  United  States  Air  Force.  In  addition 
to  the  standard  meteorological  instrumentation,  this  station  also  has  a  cloud 
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height  radar,  AN/TPQ-li,  operating  routinely.  The  fourth  requirement  was 
met  by  the  inclusion  of  standard  narrow-  and  wide-angle  pyrheliometei  s*  in 
the  experimental  setup. 

B.  INSTRUMENTATION 

The  basic  components  ol  the  instrument  designed  to  measure  the 

angular  scattering  profiles  of  clouds,  hereafter  called  the  cloud  scattering 

instrument,  are  a  lens,  a  series  of  field  stops,  an  integrating  sphere,  and 

a  detector.  The  arrangement  is  shown  in  Figure  1.  The  lens  is  a  Nikkor 

21mm  local  length,  f/4  35mm  camera  lens.  The  field  stops  are  thirteen 

holes  drilled  in  a  rotating  disc  plus  an  opaque  position.  The  hole  sizes 

correspond  to  field  of  view  ranging  from  0.  95°  to  32.  2°  half  angle.  The 

integrating  sphere  is  made  from  two  9-inch  diameter  plastic  hemispheres 

with  holes  drilled  at  their  "poles"  for  the  entering  and  exiting  light.  The 

sphere  is  painted  on  the  inside  with  3M  "White  Velvet  Coating",  which  is  a 

50-50  mixture  of  titanium  dioxide  and  silicon  dioxide.  The  detector  is 

prevented  from  seeing  direct  radiation  by  a  spider-mounted  baffle  on  the 

optic  axis  1.  5  inches  from  Ihe  face  of  the  detector.  The  detector  is  an 

RCA  7764  multiplier  phototube  with  an  S-ll  response.  This  tube  has  its 

C  °  . 

maximum  response  between  3900A  and  4900A,  and  50%  sensitivity  points 
at  3500  and  5700A  (Ref.  11). 

The  instrument  is  enclosed  in  a  weather-tight  box  with  an  outer  window 
of  glass  with  a  thin  layer  of  evaporated  chromium  to  act  as  a  neutral  density 
filter  and  heat  reflector.  Mounted  on  the  same  base  plate  and  bore  sighted 


These  instruments  were  loaned  to  The  Aerospace  Corporation  by  the  Army 
Electronics  Proving  Ground  at  Fort  Huachuca,  Arizona,  through  the  efforts 
of  Capt.  L.  Johnson  of  Detachment  50,  6th  Weather  Wing,  Space  and  Missile 
Systems  Organization,  Air  Force  Systems  Command,  Los  Angeles,  California 
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with  the  cloud  scattering  instrument  is  an  Epply  narrow-angle  (2.85°  half 
angle  field  of  view),  normal  incidence  pyrheliometer .  The  pyrheliometer 
has  a  visual  indicator  of  correct  solar  pointing  that  can  be  used  when  the  sun 
is  shining  to  adjust  the  alignment  of  the  two  instruments  during  operation. 

The  instrument  and  the  narrow-angle  pyrheliometer  are  mounted 
on  a  conventional  astronomical  telescope  clock  drive  whose  polar  axis 
was  carefully  aligned  to  true  north.  An  appropriate  declination  adjustment 
of  the  instruments  permitted  correct  tracking  of  the  sun  with  the  sun 
centered  in  the  fields  ot  view  of  the  apertures.  Electrical  and  signal 
connections  are  made  through  a  set  of  slip  rings  mounted  on  the  polar 
axle  to  prevent  wind  up  of  the  cabling  by  the  rotation  of  the  drive. 

The  instrument  was  placed  on  a  wooden  platform  on  the  roof  of  a 
hangar  at  Eglin  Air  Force  Base  approximately  50  feet  above  the  ground. 
Figures  2  and  3  show  the  instrument  in  place.  An  Eppley  wide-angle 
(2 it  field  of  view),  horizontal  incidence  pyrheliometer  was  placed  in  an 
unshaded  location  on  the  roof  of  a  small  penthouse  near  the  instrument. 

The  power  and  signal  cabling  was  run  approximately  200  feet  to 
the  Weather  Operations  room  on  the  floor  of  the  hangar  where  the  record¬ 
ing  electronics  were  housed.  Because  of  the  long  cable  run,  preamplifiers 
were  included  in  the  weather-proof  box. 

The  recording  electronics  consisted  of  a  Hewlett-Packard  Model 
3440A  Digital  Voltmeter  with  a  Model  3443A  High  Gain/Auto  Range  plug-in 
unit,  a  Hewlett-Packard  Model  562AR-J74  Digital  Recorder,  a  Parabam 
Digital  Clock/ Calendar,  and  additional  automatic  sequencing  electronics 
built  at  The  Aerospace  Corporation. 

The  operation  of  the  experiment  was  as  follows.  A  24-hour  timer 
turned  the  electronics  on  before  sunrise  and  off  after  sunset.  A  set  of  data 
8  taken  automatically  every  10  min  during  this  time.  The  instrument 
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Figure  2.  Photograph  o£  the  Instrument  as  Seen  From  the  Ground 
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Photograph  of  Instrument  as  Seen  From  the  Roof 
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sequenced  through  the  thirteen  apertures  from  smallest  to  largest,  stopping 
1-2  seconds  at  each  aperture  to  allow  recording  of  the  signal.  After  sequenc¬ 
ing  through  the  thirteen  apertures,  the  instrument  turned  itself  off  at  the 
opaque  position.  A  zero  (opaque)  reading  was  taken  at  the  start  and  end  of 
each  sequence.  After  cycling  through  the  apertures,  the  instrument 
recorded  the  signals  from  the  two  pyrheliometers .  Thermistors  located  in  the 
cloud  instrument  and  each  of  the  pyrheliometers  were  also  sampled.  Finally, 
the  minute,  hour,  and  day  number  v/ere  recorded  from  the  digital  clock/ 
calendar.  The  entire  cycle  took  about  a  minute,  after  which  the  electronics 
went  into  standby  until  ten  minutes  elapsed. 

The  pyrheliometric  measurements  were  made  with  two  standard 
Eppley  pyrheliometers,  a  50 -junction  hemispheric  pyrheliometer,  and  a 
Model  15  normal -incidence  pyrheliometer  (Ref.  12).  The  cloud  height  radar 
in  operation  at  Eglin  Ai~  Force  Base  was  a  Model  AN/TPQ-11  Radar  Cloud 
Detecting  Set  operating  in  the  frequency  band  (33-36  GHz)  (Ref.  13).  The 
facsimile  chart  recordings  from  this  instrument  were  provided  to  The 
Aerospace  Corporation  personnel  who  made  the  interpretations.  These 
interpretations  were  made  by  combining  the  standard  observations  of  the 
clouds  and  the  indicated  clouds  on  the  charts.  This  was  necessitated  by 
the  inability  of  the  AN/TPQ-11  to  detect  all  clouds.  Kantor  (Ref.  14)  has 
analyzed  the  detection  capabilities  of  the  AN/TPQ-11  radar  and  found  that 
only  62%  of  visually-observed  overcast  clouds  were  detected  by  the  radar. 

Standard  meteorological  observational  data  were  supplied  by  the 
personnel  at  Eglin  Air  Force  Base  on  WBAN  Form  10  "Surface  Weather 
Obse  rvations.  "  The  WBAN  Form  10  includes  the  following  data  used  in 
this  experiment:  dme,  sky  condition,  horizonta  visibility,  temperature 
and  dew  point,  total  sky  cover,  remarks,  and  supplemental  coded  data. 
Observations  were  made  hourly  except  when  rapidly  changing  conditions 
warranted  more  frequent  observations.  The  supplemental  coded  data 
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included  the  synoptic  code  observations  of  low,  middle,  and  high  clouds 

(1C,  C.  C  )  (Ref.  15).  These  were  recorded  every  three  hours.  These  synoptic 
Li  M  H 

measurements  were  used  to  identify  the  clouds  present  during  our  mea¬ 
surements.  Identifications  between  the  synoptic  observations  were  made 
by  interpolating  between  the  bracketing  observations  using  the  remarks 
recorded  by  the  observer  as  to  advancing  or  receding  clouds. 

C.  CALIBRATION 

The  calibration  of  the  cloud  scattering  instrument  entailed  several 
steps.  All  the  measurements  were  made  in  a  darkroom  using  a  small  but 
very  bright  tungsten  coil  lamp  mounted  in  the  same  horizontal  plane  as  the 
cloud  scattering  instrument  and  about  twenty  teet  away,  A  slight  adjust¬ 
ment  in  the  focus  ot  the  lens  of  the  instrument  was  made  to  correct  for 
the  non  parallelism  of  the  light  from  twenty  feet. 

The  normalization  factors  between  the  apertures  were  determined 
by  measuring  the  on-axis  signal  for  each  aperture  and  normalizing  to  the 
smallest  aperture  (1°).  These  normalization  factors  are  listed  in  Table  !• 

The  effective  angular  field  of  view  of  each  aperture  was  determined  as 

4 

follows.  With  no  limiting  aperture,  the  intensity  is  given  by  a  cos  0  angular 
2 

distribution.  The  signal  was  measured  at  a  number  of  angles  on  both 

sides  of  the  optic  axis  (in  a  horizontal  plane  only).  The  two  points  at  which 

4 

the  signal  dropped  to  50%  of  the  above -mentioned  cos  0  value  were  averaged 
to  give  the  effective  field  of  view  of  a  given  aperture.  The  results  of  these 
measurements  are  given  in  Table  2.  The  areas  under  these  curves  of 
angular  response  were  measured  with  a  planimeter  to  determine  the  solid 
angle  factor  for  each  aperture.  These  solid  angle  factors  are  given  in 
Table  3. 


The  cos  6  represents  the  fall  off  of  intensity  of  off-axis  rays  due  to 
geometric  factors.  See,  for  example:  Jenkins  and  White,  Fundamentals 
of  Optic 8,  Chapter  7. 
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Table  1.  Apei  .ure  Normalization  Factors 


Aperture 

Normalization  Factor 

1 

(si.  0000) 

2.  5 

0. 9882 

5 

0. 9824 

7 

0. 9824 

10 

0. 9719 

12 

0. 9644 

14 

0. 9516 

16 

0.  9143 

18 

0.  8947 

22 

0. 8330 

26 

0. 8029 

30 

0. 7928 

32 

0. 7966 
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Table  2.  Effective  Fields  of  View 


Aperture 

Effective  Field  of  View 
(Half  Angle) 

1 

0.  9583° 

2.  5 

2. 400° 

5 

4.  858° 

7 

7. 171° 

10 

9. 763° 

12 

12. 076° 

14 

14.  290° 

16 

16. 198° 

18 

18. 156° 

22 

21. 882° 

26 

25.  882° 

30 

30.042° 

32 

32. 170° 

•i 


Table  3.  Solid  Angle  Factors 


Aperture 


Solid  Angle  Factor 


1 

2.  5 
5 
7 

10 

12 

14 

16 

18 

22 

26 

30 

32 


0. 0002688 
0. 0017617 
0. 0071867 
0. 015631 
0.028955 
0. 044267 
0. 061882 
0. 079400 
0. 099587 
0. 14410 
0. 20062 
0. 26868 
0. 30705 


3 

The  conversion-calibration  factor  in  langleys/min/volt  for  each  of 
the  pyrheliometers  is  given  by  the  Eppley  Company  with  their  instruments. 

An  additional  calibration  factor  required  is  the  gain  of  the  preamplifiers 
that  take  the  low-level  pyrheliometer  output  and  amplify  it  prior  to  its 
being  transmitted  down  the  long  signal  cable.  It  was  impossible  to  cali¬ 
brate  the  entire  pyrheliometer-preamplifier-signal  cable  assembly  as 
installed  in  Florida.  An  indirect  method  of  calibration  had  to  be  applied. 

Several  crystal-clear  days  with  very  high  reported  horizontal 
visibility  occurred  in  the  early  period  of  the  experiment.  The  calculated 
transmittance  of  a  Rayleigh  atmosphere  in  the  spectral  region  passed  by 
the  glass  bulb  of  the  pyrheliometers  is  0.  904  and  the  calculated  trans¬ 
mittance  of  a  turbid  atmosphere  is  0.718.^  Since  the  atmospheric  con¬ 
ditions  on  these  crystal-clear  days  seemed  to  be  clearer  than  those 
described  by  Elterman  for  his  "turbid"  atmosphere,  the  simple  Rayleigh 
atmospheric  transmittance  of  0.  90  was  adopted  for  the  pyrheliometric 
readings  on  these  days.  With  the  use  of  this  value  of  the  transmittance,  the 
calibration  factors  of  the  pyrheliometers  were  calculated.  The  absolute  accuracy 
of  the  pyrheliometric  measurements  does  not  affect  the  shapes  of  the 
angular  scattering  profiles  that  are  more  significant  to  this  study  than 
the  absolute  magnitudes. 


A  langley  is  1.  0  cal/cm^. 

^The  calculated  transmissions  are  based  on  Elterman's  atmospheric 
attenuation  model  as  given  in  Handbook  of  Geophysics  and  Space  Environ¬ 
ments,  U.  S.  Air  Force  Cambridge  Research  Laboratories,  hlew  York: 
McGraw-Hill  Book  Co  ,  Inc.  ,  1965,  Chapter  7.  The  spectral  range  was 
considered  to  be  that  passed  by  the  glass  window  or  envelope  of  the 
pyrheliometers;  namely  0.  35  to  2.  5  microns.  The  sun  with  a  spectral 
distribution  as  given  in  Chapter  16  of  the  Handbook  was  used  for  the 
radiation  source. 
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IV.  EXPERIMENTAL  OPERATION  AND  DATA  SAMPLE 


A.  PERIOD  OF  EXPERIMENT 

The  experiment  was  operated  between  9  February  and  18  October 
1968.  Because  of  equipment  malfunctions,  data  were  obtained  only  for  the 
periods  9  February  through  3  April,  5  April  through  9  April,  17  April 
through  22  May,  12  June  through  21  July,  28  August  through  4  October, 
and  16  October  through  18  October.  The  data  for  15  May  through  22  May 
were  lost  in  the  mails.  The  narrow-angle  pyrheliometer  was  inoperative 
for  the  period  12  June  through  21  July.  The  data  for  the  period  28  August 
to  the  end  of  the  experiment  are  considered  unreliable  because  of  a 
belatedly  discovered  possible  misfocus  and  f/stop  change.  Reliable, 
complete  data  are  available,  therefore,  only  for  88  days  between 
9  February  and  14  May.  hi  addition,  incomplete  data  (narrow-angle 
pyrheliometer  inoperative)  are  available  for  the  40-day  period  between 
12  June  and  21  July. 

B .  DATA  SAMPLE 

During  the  88-day  period  for  which  complete  data  are  available,  a 
total  of  3296  measurements  were  made  during  times  when  the  sun  was 
higher  than  35°  above  the  horizon.  Of  these  3296  measurements,  530 
were  made  during  completely  clear  skies  (0/10's  cloud  cover)  and  1200 
during  completely  overcast  skies  (10/10's  cloud  cover).  The  1566  mea¬ 
surements  made  durirg  periods  of  broken  or  scattered  clouds  were  not 
used  for  further  analysis  because  of  the  uncertain  location  of  holes  in  the 
clouds  with  respect  to  the  sun's  position.  In  addition,  measurements  of 
overcast  clouds  made  when  more  than  one  layer  of  clouds  was  reported 
were  also  not  used  for  further  analysis.  This  war  because  of  uncertainties 
about  breaks  in  the  individual  layers.  This  left  a  sample  of  530  measure¬ 
ments  during  completely  clear  skies  and  398  measurements  made  during 
single-layer,  completely  overcast  skies.  Of  the  398  overcast  measure¬ 
ments,  64  were  discarded  because  of  various  problems  with  the  data. 
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Table  4  gives  a  breakdown  of  the  32  96  measurements  by  che  amount 
of  cloud  cover.  Table  5  gives  a  breakdown  of  the  1200  overcast  measurements 
by  the  number  of  cloud  layers  reported.  Table  6  further  breaks  down  the 
single-layer  overcast  measurements  into  reported  cloud  types.  This  table 
gives  the  breakdown  of  the  measurements  upon  which  the  an.ilysis  was  based. 


Table  4.  Number  of  Measurements  Made  by  Amount  of  Cloud  Cover 

(Solar  Elevation  Anglo  2  35°) 


Cloud  Cover 

Number  of  Measurements 

_  i 

0/10 

530 

1/10 

236 

2/10 

152 

3/10 

142 

4/10 

189 

5/10 

109 

6/10 

112 

7/10 

215 

8/10 

214 

9/10 

168 

10/10 

1200 

Broken  (6/10  to  9/10) 

7 

Scattered  (1/10  to  5/10) 

16 

Unknown 

6 

Total 

1 

3296 
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Table  5.  Number  of  Measurements  of  10/10's  Cloud  Cover  by 
Number  of  Cloud  Layers  Reported 


Number  of  Cloud  Layers 

Number  of  Measurements 

1 

398 

2 

486 

3 

249 

>  4 

67 

Total 

1200 
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Table  6.  Number  of  Complete  Measurements  of  Single -Layer 
10/10's  Cloud  Cover  by  Cloud  Type 


Cloud  Type 

Name 

No.  of 

Measurements 

Low  clouds 

CI  =3 

Cumulonimbus 

3 

O 

F 

n 

4^ 

Stratocumulus  cumulogenitus 

5 

CL  =  5 

Stratocumulus 

34 

Cl  =  6 

Stratus 

127 

Middle  clouds 

CM=  1 

Altostratus 

5 

CM  ~  3 

Altocumulus  translucidus 

5 

CM  =  7 

Altocumulus  opacus 

17 

High  clouds 

CH=  1 

Cirrus  fibratus 

5 

CH=7 

Cirrostratus 

133 

Incomplete  data 

64 

or  unknown  cloud 
type 

V.  DATA  REDUCTION 


The  first  stage  in  the  data  reduction  scheme  was  the  laborious  manual 
transcription  of  the  data  from  the.  digital  printer  tapes  to  computer  coding 
sheets.  This  was  done  by  the  computress  staff  in  the  Space  Physics  Laboratory 
at  The  Aerospace  Corporation. 

Three  cards  were  vised  to  encode  all  the  data  for  each  observation. 

The  data  included  for  each  observation  were  the  following:  day,  time,  wide- 
angle  pyrheliometer  reading,  narrow-angle  pyrheliometer  reading,  photo¬ 
tube  initial  and  final  zero  readings,  the  thirteen  aperture  readings,  the 
thermistor  readings  of  the  two  pyrheliometers  and  the  phototube,  an 
abbreviated  description  of  the  weather  and  obscurations  to  vision,  the 
total  cloud  cover,  the  horizontal  visibility,  the  air  temperature  and  dew 
point,  and  a  description  of  up  to  four  layers  of  clouds.  This  description 
included  the  coverage  of  the  individual  layer,  the  type  of  cloud,  and  the 
base  and  top  of  the  layer. 

A  computer  program  was  written  (by  R.  T.  H.  )  which  took  these 
three  cards  for  each  observation  and  computed  the  mean  relative  radiance 
and  irradiance,  the  wide-  and  narrow-angle  transmittances,  the  relative 
humidity,  and  the  solar  elevation  angle.  All  these  data  were  printed  by 
the  computer  on  one  page  together  with  verbal  descriptions  of  the  weather 
and  clouds.  Plots  were  made  on  3  5-mm  film  of  the  mean  relative  radiance 
and  irradiance  as  a  function  of  the  field  of  view.  An  example  of  a  page  of 
computed  data  is  given  in  Figure  4. 
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The  relationships  used  to  compute  the  various  output  quantities 


are  the  following: 

IRRAD(I)  =  [READING^)  -  ZERO]  x  NFAC(I)  (1) 

RAD(I)  =  IRRAD{I)/SAF(I)  (2) 

RELHUM  =  [VP(DEWPT)/VP(TEMP)]  x  100  (3) 

sin  V  =  sin  6  sin  4>  +  cos  6  cos  ft  cos  $  (Ref.  16)  (4) 

ft  =  TIME  -  11.76722  (5) 

cos  P  =  sin  V  [cos  e  +  sin  t  cos  A  tan  4>]  -  (Ref.  16)  (6) 

sin  6  sin  e  cos  A  sec  <j>  +  cos  6  sin'  ft  sin  t  sin  A 
I  =  2.000/Rv2  (Ref.  16)  (7) 

TwW  =  ^WA/U 

Tw  =  TwW/cOS  ^  (Ref‘  16)  {9) 

Tn  =  [NA/lo]  (10) 


The  definitions  of  the  quantities  in  these  equations  are: 

IRRAD(I)  =  mean  relative  irradiance  for  aperture  I 
ZERO  =  average  of  initial  and  final  zero  readings 
READING(I)  =  reading  for  aperture  I 
NFAC(I)  =  normalization  factor  for  aperture  I 
RAD(I)  =  mean  relative  radiance  for  aperture  I 
SAF(I)  =  solid  angle  factor  for  aperture  I 
RELHUM  =  relative  humidity 

VP(I)  =  vapor  pressure  of  water  at  temperature  I 
DEWPT  =  dew  point  temperature 
'5 

Not  all  these  quantities  were  actually  used  in  the  final  analysis. 
^Printouts  of  the  complete  output  data  a  *e  available  from  the  author. 

It  should  be  noted  that  cos  P  =  sin  y  for  a  horizontal  plane,  and 

Equation  9  would  take  its  usual  form  T  =T  ,,,/sin  v 

w  wW'  Y  - 
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TEMP  =  air  ‘emperature 

Y  =  solar  elevation  angle 

6  =  solar  declination  angle 

9  ~  latitude  (30°29'N  for  Eglin  AFB) 

£2  =  hour  angle  (11.  76722  is  the  local  Central  Standard  Time  of 
true  solar  noon  at  Eglin  AFB) 

P=  angle  between  the  solar  elevation  angle,  Y  .  and  the  normal 
to  an  arbitrarily  inclined  plane 

<  =  angle  between  a  horizontal  plane  and  an  arbitrarily  inclined 
plane 

A  =  azimuth  of  the  vertical  plane  containing  the  normal  to  an 
arbitrarily  inclined  plane 

Iq  =  amount  of  solar  radiation  incident  at  the  top  of  the  atmosphere 
(2.000  langleys  per  minute  is  the  solar  constant  used) 

=  radius  vector  of  the  sun  (earth-sun  distance) 

TwW  =  wide -angle  transmittance  of  the  sun's  radiation  not  taking 
into  account  the  solar  elevation  angle 

T  =  narrow-angle  transmittance 
w  ° 

WA  =  wide-angle  pyrheliometer  calibrated  reading 

NA  =  narrow-angle  pyrheliometer  calibrated  reading 

The  quantities  (3  and  t  above  were  necessitated  by  a  small  tilt 
(~3°WNW)  of  the  wide-angle,  horizontal  incidence  pyrheliometer  which 
occurred  shortly  after  its  installation.  The  correction  resulted  in 
approximately  a  5%  change  in  the  wide-angle  pyrheliometer  readings. 
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VI.  ANALYSIS  AND  RESULTS 


A.  ANALYSIS 

The  goal  of  the  experiment  was  the  development  of  a  realistic 
model  of  the  angular  scattering  profiles  of  clouds  for  use  in  designing 
visible-light  optical  systems.  It  was  hoped  that  the  observed  scatter¬ 
ing  profiles  could  be  correlated  with  some  routinely  observed  meteor¬ 
ological  quantity  or  quantities  so  that  system  performances  could  be 
predicted  from  readily  available  meteorological  data. 

Various  correlations  and  relationships  were  tiled.  The  most 
obvious  is  a  correlation  between  the  irradiance  measurements  and  the 
cloud  type.  This  correlation  was  tried  and  failed.  This  failure  is 
probably  due  to  the  great  variability  in  cloud  characteristics  within  a 
cloud  type.  The  same  cloud  can  exist  over  a  wide  range  of  thicknesses 
and  opacities. 

The  failure  of  this  obvious  correlation  prompted  a  search  for  some 
other  way  of  describing  the  scattering  profiles,  such  as  a  combination 
of  the  measures  of  the  direct  and  diffuse  solar  radiation.  This  descrip¬ 
tion  should  be  reasonable  in  light  of  the  observation  that  on  a  clear  day 
the  sun's  radiation  is  virtually  entirely  direct  and  that  on  a  heavily 
overcast  day  the  sun's  radiation  which  penetrates  to  the  earth  is  entirely 
diffuse.  All  other  cases  should  be  some  combination  of  these  two  extremes. 
The  intermediate  cases  can  be  described  by  a  bright  spot  at  the  position 
of  the  sun  surrounded  by  a  diffusely  scattering  region. 

The  recults  as  derived  by  the  computer  analysis  are  in  terms  of 
relative  irradiances.  A  more  useful  quantity  is  the  transmittance.  Since 
there  was  no  wav  to  derive  absolute  transmittances  from  the  relative 
irradiances  an  alternative  was  to  derive  the  transmittances  relative  to  a 
clear  sky.  The  transmittances  could  be  determined  from  the  relative 


-29- 


Figure  5.  Relative  Irradiance  Plot  for  1205  Hours,  1  March  1968 


Table  7.  Relative  Irradiance  and  Least-Squares  Values  for  Clear  Sky  Case 


e 

Relative  Irradiance 

Least-Squares  Value 

0. 9583 

— 

2.06447 

2.  400 

2.  06147 

2. 06509 

4.  858 

2. 06995 

2. 06613 

7.  171 

2. 06613 

2. 06712 

9.  763 

2. 06845 

2. 06822 

12. 076 

2.  06969 

2. 06920 

14. 290 

2.  07364 

2. 07014 

16.  198 

2. 07244 

2. 07096 

18. 156 

2.06909 

2. 07179 

21.  882 

2. 07186 

2. 07337 

25.  882 

2.07277 

2. 07508 

30.  042 

2.07434 

2.07685 

32. 170 

2.08186 

2. 07775 

(Relative  Irradiance)  .  .  =  2.064064  +  0.00042550 

'clear  sky 

F;gure  6a.  Sample  Plot  Showing  Example  of  Quality  of  Fit 


Figure  6c.  Sample  Plot  Showing  Example  of  Quality  of  Fit 


Figure  6d.  Sample  Plot  Showing  Example  of  Quality  of  Fit 


0.40 


<x> 


3DNVlllWSNVdx 


-37- 


Figure  6e .  Sample  Plot  Showing  Example  of  Quality  of  Fit 


With  the  coefficients  A  and  B  available  as  descriptors  of  the 

transmittance  profiles,  correlations  of  these  coefficients  with 

routinely  measured  actinometric  quantities  could  be  tested.  The 

four  combinations  of  A  and  B  with  T  and  T  were  tested  for  possi- 

w  n  r 

ble  correlations.  Separate  tests  were  made  for  "good",  "fair", 

"poor",  "good  +  fair",  and  all  data.  Some  of  these  plots  are  shown 
in  Figures  7-10. 

Looking  at  Figure  7  there  appears  to  be  a  linear  correlation 
with  an  origin  intercept  between  and  A.  This  is  reasonable  in  view 
of  the  fact  that  both  are  measures  of  the  direct  sunlight.  A  linear 
correlation  with  origin  intercept  between  T^  and  B  is  again  indicated 
in  Figure  8  although  the  scatter  is  fairly  great  at  large  values  of  T  . 
This  correlation  is  also  reasonable  in  view  of  the  fact  that  both  quanti¬ 
ties  are  measures  of  the  diffuse  radiation. 

Figure  9  shows  no  simple  correlation  between  Tfl  and  B  as 

would  be  expected.  Figure  10  illustrates  a  more  complex  situation.  A 

linear  correlation  seems  to  exist  between  T  and  A  but  with  a  non-zero 

w 

intercept.  It  would  seem  that  a  threshold  value  of  T  must  be  exceeded 
r  w 

before  A  has  a  nonzero  value.  This  is  reasonable  from  the  observation 
that  the  sun  can  be  seen  through  the  clouds  only  for  a  range  of  transpar¬ 
ency  of  the  clouds.  For  more  opaque  clouds  no  direct  sunlight  can  be 
seen. 

Linear  least-squares  fits  of  these  data  and  their  correlation 
coefficients  were  calculated  to  quantify  these  visual  correlations.  The 
"poor"  data  were  omitted  from  these  calculations.  The  calculations  for 
Tw  versus  B  and  Tn  versus  A  are  straightforward.  The  calculation  of 
Tw  versus  A  requires  that  the  threshold  point  be  determined  first.  The 
threshold  point  was  found  by  using  only  these  points  in  the  least-squares 


-38- 


calculations  whose  coordinate  was  above  some  arbitrary  value. 

The  threshold  point  was  defined  when  the  cutoff  most  nearly 
equalled  the  T^-axis  intercept.  The  threshold  determined  by  this 
method  was  0.  54.  The  results  of  the  linear  least-squares  calcula¬ 
tions  are  plotted  in  Figures  11-13  and  tabulated  in  Table  8. 

An  additional  correlation  between  T  and  T  was  tested. 

w  n 

These  data  are  shown  in  Figure  14.  The  results  are  qualitatively 

similar  to  those  for  the  correlation  between  T  and  A.  The  least- 

w 

squares  fit  with  a  Tw  threshold  of  0.  54  gave  the  straight  line  shown 
in  Figure  14  and  the  results  tabulated  in  Table  8. 

B.  RESULTS 

In  light  of  the  good  correlation  between  T  and  A  and  between 
Tw  and  B,  it  would  appear  that  a  reasonable  model  of  the  angular 
scattering  characteristics  of  clouds  has  been  found  which  can  be  quan¬ 
titatively  calculated  from  two  routine  meteorological  observations. 
Combining  the  results  tabulated  in  Table  8  with  Equation  11  results 
in  the  following  formula  for  the  visible  light  angular  transmittance  of  a 
cloud  relative  to  that  of  the  clear  sky  as  a  function  of  the  narrow-  and 
wide-angle  pyrheliometric  transmittance: 

T  =  -0.  0056  +  0.  89  X  T  +  [-0.028  +  1.09  X  T  J  X  sin20  (12) 

n  *-  w 

The  formula  is  valid  at  least  for  0  ^  35°.  For  largei  values  of  0  an 
extrapolation  of  unknown  accuracy  is  assumed. 

How  well  this  model  will  predict  the  angular  transmittance  pro¬ 
files  is  shown  in  Figures  15a-e.  The  data  are  the  same  as  used  in 
Figures  6a-e.  The  lines  in  the  figures  are  calculated  from  Equation  12; 
the  points  are  the  observed  data.  The  agreement  of  the  calculated  lines 
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Figure  12.  Plot  of  T  vs  A  Showing  Linear  Least-Squares  Fitted  Line 


Figure  13.  Plot  of  Tw  \  s  B  Showing  Linear  Least-Squares  Fitted  Line 
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Figure  15a.  Sample  Plot  Showing  Accuracy  of  Model 


Showing  Accuracy  of  Model 


igure  1  5d.  Sample  Plot  Showing  Accuracy  of  Model 


Figure  15e.  Sample  Plot  Showing  Accuracy  of  Model 


with  the  observed  points  is  quite  satisfactory.  An  estimate  of  the 
accuracy  of  a  calculated  transmittance  is  ±0.02  or  ±10%>  whichever 
is  greater. 


VII.  INTERPRETATION  AND  APPLICATION 


A.  INTER]  .*  STATION 

The  success  of  the  anguljir  scattering  model  developed  here 
is  really  a  validation  of  the  idea  that  the  transmission  of  light  through 
the  atmosphere  has  two  components:  a  direct  component  and  a 
scattered  component,  and  that  these  two  components  are  separable. 

The  relative  magnitudes  of  the  two  components  varies  with  the  character¬ 
istics  of  the  scattering  media  in  ^he  light  beam.  The  di.ect  component 
is  important  only  when  the  wide  angle  transmittance  of  the  scattering 
media  is  greater  than  roughly  0.  50. 

Diffuse  radiation  is  usually  taken  to  mean  an  intensity  decrease 
with  the  cosine  of  the  angle  to  the  normal  of  the  surface  from  which  the 
radiation  emerges.  Such  a  surface  appears  equally  bright  from  all 
directions  because  this  cosine  dependence  just  compensates  for  the  fore¬ 
shortening  of  the  surface  as  seen  by  the  observer.  The  diffuse  component 
2 

of  T,  B  sin  6,  which  was  calculated  from  this  experiment  would  have 
been  simpler  to  explain  if  the  measurements  could  have  been  made  by 
observing  at  the  zenith  so  that  the  vector  to  the  sun  would  have  been 
perpendicular  to  the  cloud  layer.  Then  G  would  have  been  the  ray  angle 
from  the  cloud  normal.  However,  the  minimum  solar  zenith  angle 
during  the  experiment  was  about  12°  and  measurements  were  necessary 
at  still  larger  zenith  angles  (to  55°)  to  obtain  a  reasonable  data  sample. 
Because  a  diffuse  surface  appears  equally  bright  in  all  directions, 
the  readings  obtained  for  the  diffuse  components  are  still  properly 
made.  Note,  however,  that  tor  nonzero  solar  zenith  angles,  the  path 
length  of  the  radiation  through  the  cloud  increases  and  the  transmittance 
decreases.  This  is  borne  out  qualitatively  by  the  diurnal  plots  of  T^ 
(denoted  PPT)  in  Ref.  1. 


It  should  be  possible  to  relate  T  and  T  directly  to  A  and  B. 
r  w  n  ' 

Indeed  T  =  T  for  8  =  90°  and  T  =  T  for  8=2.  85°.  If  we  use  these 
w  n 

relations  to  evaluate  A  and  B  we  find  that: 

T  »  T  +  (T-  -  T  )  X  sin28 
n  '  w  n' 

These  values  of  A  and  B  do  not  correlate  with  the  experimental  data  as 
well  as  those  deduced  above.  Possible  partial  explanations  for  this 
discrepancy  include: 

1.  The  cloud  scattering  instrument  measures  in  the  spectral 
band  from  abouf  0.  35  to  0.  57  microns  with  an  S-ll  photomultiplier 
respor.se.  The  pyrheliometers  are  thermal  detectors  which  respond  to 
all  wavelengths  transmitted  by  the  atmosphere  and  by  their  glass  envel¬ 
opes  or  windows,  i.  e.  from  about  0.  35  to  2.  5  microns.  The  scattering 
properties  of  the  cloud  droplets  are  wavelength  dependent. 

2.  The  corrections  for  atmospheric  absorption  used  in  the 
normalizations  are  not  the  same  for  the  two  types  of  instruments. 

3.  The  model  was  derived  only  for  8  <  35°.  An  extrapolation 

of  the  model  to  8  =  90°  is  of  unknown  accuracy. 

B.  APPLICATION 

The  model  derived  for  the  angular  scattering  profile  of  visible 
solar  radiation  through  clouds  has  several  possible  limitations  on  its 
validity  and  applicability.  The  primary  limitation  is  that  the  model  is 
valid  only  for  single,  overcast  cloud  layers.  Multiple  cloud  layers 
cannot  be  treated  for  the  simple  reason  that  the  light  impinging  on  the 
second  layer  after  emerging  from  the  first  layer  no  longer  can  be  con¬ 
sidered  to  be  collimated.  Rather  it  has  a  scattered  component  that  will 
diffusely  illuminate  the  second  layer,  thereby  grossly  deviating  from  the 
conditions  of  the  model. 
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Likewise  the  model  is  inapplicable  for  broken  clouds  even 
with  the  source  behind  one  of  the  clouds.  The  direct  component  is 
affected  by  the  presence  of  broken  clouds  to  the  same  extent  as  it  was 
for  an  overcant  laye? ,  However,  the  scattered  component  will  be 
grossly  changed  because  of  the  non-uniform  scattering  medium.  The 
change  can  be  either  to  increase  or  decrease  the  apparent  scattered 
component.  The  i.ncrease  could  result  from  strong  reflections  off  the 
sides  of  the  clouds  into  the  field  of  view  of  the  sensor.  A  decrease 
might  result  from  a  fortuitous  arrangement  of  the  clouds  sc  as  to 
preferentially  scatter  the  beam  away  from  the  sensor  with  little  or  no 
scattering  back  into  the  beam.  Such  a  case  might  be  a  single  small 
cloud  directly  in  front  of  the  source  with  no  other  clouds  in  the  sky. 

A  further  limitation  is  that  the  model  is  applicable  only  to 
collimated  light.  Diverging  or  converging  light  might  deviate  wildly 
from  the  results  predicted  by  the  model.  This  aspect  of  the  model  is 
subject  to  experimental  verification  albeit  with  considerable  difficulty. 

What  has  been  developed  here  is  a  model  which  describes  the 
angular  scattering  characteristics  of  collimated  visible  light  through 
a  single,  overcast  layer  of  clouds.  The  model  gives  satisfactory 
agreement  with  experimental  observations.  The  model  should  be  of 
great  use  in  the  evaluation  and  prediction  of  the  performance  of  visible 
light  optical  systems  in  the  presence  of  clouds.  It  should  also  serve 
as  a  useful  stepping  off  point  for  the  development  of  models  of  more 
complex  situations  such  as  broken  clouds  or  multiple  cloud  layers.  It 
should  also  serve  as  a  useful  test  for  the  theoretical  models  of  collimated 
light  incident  on  an  infinite  scattering  layer. 
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APPENDIX 


EGLIN  DATA  FOR  10/10’s  COVER  SORTED  BY  CLOUD  TYPE 


A  listing  of  all  334  data  samples  is  shown  on  the  following  pages. 
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eGLIN  CLOD  DATA  RP  10/10THS  CO V£R  SORTED  BY  CLOUD  TYPE 
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r*:»j3<co^crx)flOot>-tvoror^4-4,^,(T'4r>vDCMroNUN4^K(VJiA»Dj'\OJ,ou>o'»Dmu>rv<r«wj*fo 
*Hin^'r><Vf\JK.(J,CJO<\jt\J.**-<*HS.,Hrr,  u>.*erv*U'\«TC»«Olf»OCir»tfK.T*Nf'>fr:S.*HtANli><V.!tX> 
|’O^J(^e^^^DU^^^'r^lfMJ^^•tfl0^’  ^OhO'rv.'CT^O.tNMPK 
tfi  *,  C  O-  in  ep  ^  (T  J  D  r  rv.  lf»  iCTePiDOf0^4,'7'C^NN«ON.HOU\in«)NJ’^(r>OMftN 


v0(MiiN®cr®v0U3^o(T'f0K®^\0f0fl0^cvj*-<a'r0fvj«0»^i0<viwot-i*Hirtcvjpwinajf^wov0j% 

(\JO'N.roa3^CViin»-i®U>OO^^^lAWflOvDfOUNm^C^OflOr^^fON.(T'CT'eOaOJ)*-»^^t€C4*fVisO 

trtrr>^-v0i/N^K^r^«Ha>ir\uN^N.^u>ii>^^KvO4>UN^irir>iAiA«Hro«£)\O<vjirvKrwUMn-^’Nm 


u>vO-tcr^oo^a<HNO'o(Vj  j'^ffrN.arf\jiTij-«'^roocrrv(T<C'crj-  -*t\jir\cvj«oj’S.r>*H«n 
OCC~(V4-«-<aO<\eOv£>r-  U\CV)K*lf»iD»-<fO»1OU><V)CrUD<r(T'^^’>o(T'W<y'0<T'K)  ^  f*>  41  K  K1  CD  S  ^  ^ 

CTlfnDf0lI>v0<rr«C^aNlAf'^(rNHOOHlf\OOCPO4'C>J<rN(\JN(T'0rrf0«0^^jraa3^fn 


aH^r0CJHf0(\JHr^Hf\jT-(a(\J^(\JfVJH(\JC'JC\JHHO^(\j»HCVJO4-OOt\J(\JC»Cj(\JDOO(V, 


(Tfl'O'ff'O'O'fl'O'O'O'oooo  o  «-*  -r-4  ^ 
SSSSN,S\SSN(eCOC0CTfl3CCa0C 


*H*-c<Hru«\i<viCVJ<'jc\JKforifofofr'  4-4,^,J‘4'ir»u\4f»iA 

«r«'cocrcDccecccCdocrcc  cCio®cr«<cffCco® 


<T'f'-irvKa3rc(vj(\j*-*<r*-»vOfNj®Lnj'N.*sX)r'>veer,*N«e‘f*r'>N.o^rof*><NJ4'inr*'r*c'<rror'>.* 

4‘^«rvr\C\£)N.«£>Kfv«Ka}«£>^-N^v£)v£>vCNN\O^Da>^>^Koa>«D^irvp>-Krr«£>rs.rs.xotnu>^\0 

of)f?nocooc'rrc>a«'eoc3rooocoorDOc'orror)roprcc'oor 

Of0^r*'r<0»DCCm«fOOfOtOHCOr4eD«^dOlMI>fH^ClAOOUM|\OOj)OOOlfl 
rl  ^  H  H  t<  r<  H  »4  rl  »•<  <H  H  H  r4  tH  r4  H  H  r<  r<  iH  H  r< 


o2C0Z0Cc0cc?Zcccc20Z07o0C'Z?oZt2Z0?7r?C00ZZ 
cKconrcrcc  o  n  c  o  o  o  e  *-o»Hr*C'oKKC'MC,MK*c,>-MK»^r“oc?>-^ 
oiooxcrciooooXloircoioioioooiXoXolXoxTixoooix 
rr  t~  j-  s  k  vr  k  »ch^Hr>-\D(D(vj»-h-^^  ®  k  k  r  k  »-  r  f- 

m  r-r*>  r  r  ru  r  e\j  rr  »o  r*  r  r  rr  k  <u  r>  ev  m  k  cv  fv* 

dd  c  r  i  c  o  o  c  r  o  •  c  o  c  o  c  •  o  #0  •  o  c  n  •  •  o  c  o  •  o  n  •  •  o  •  o  c  o  *  • 

1  c  r  o  z  r  zocrozzooczo 
ryoycvcooyyopoye 
®oa>OZJ1ooo«OXo®CT»-'J>Z\OZoX-y^®^Z^or^Xca«ZXc5X^i«r«zy 
Hn\iK)Dwr^Mr'^«3w^w^NDryDf03(\(r;rODfvjMr4jN^DOro«HwD3 

r*  c  kit  O'  k  «c  cr  •ff’cr'CCK  ^  w  it  c  cr  r  k*  C‘ff,r>«MTff  *>  <r  <r  k  k  «  k  w  f«  a 

»-<(\JK)^HfUrJOOlf\K  OO^OriO'f.  ^  r  tf%C<\J*4Aj|ftfcO*9*4 

cir^Jcorjrr  .*  c*  it  cr  j-  a  e>  cr  c*»  o(T  00O  O'o^cchoCu  «h  «h  0*  o  o  «-<  <t 


«v»  a  cr  (r 

^  O  fv.  o 


<r  r  ®c«r«*rcr«CK  <r  w  cr  c  <r  r  k-  (rtr  r>  »  nit  r  »o  o>  cr  k  k  <>  k  r>  rn  < 

r<Nrjo  otf\  k  oo^or^^r. 

c'OKj>r>K  ^c«if<rj-oe»<r^o»^o(rooa(rotfNc<wo«(y  *h  «h  O'  o  o  *h  < 

w* 

^O'crocxc\j<r«#cr<\.a'ocr  ^  cr  ^  <r  .*  <r  <r  r  *#  .*  (r  <r  ■*-«-*  «v  ru  •*,*<\jf'uC’fu<vro 

rjro^oioNc^cHC'^r  rv  o  <v  r  •*,*,*.*,*.* 

•4  r1  r4  H  v*  r4  r1  v4  H  r4  r<  H  H  H  H  »4  H  r4  v4 
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